Parity Symmetry in QED3 
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Schwinger-Dyson equations are used to study spontaneous chiral and parity symmetry breaking 
of three dimensional quantum electrodynamics with two-component fermions. This theory admits 
a topological photon mass that explicitly breaks parity symmetry and generates a fermion mass. 
We show that the pattern of symmetry breaking maintains parity but breaks chiral symmetry. We 
also find that chiral symmetry is restored at a critical number of fermion flavours in our truncation 
scheme. The Coleman-Hill theorem is used to demonstrate that the results are reasonably accurate. 

PACS numbers: 11.30.Qc, ll.15.Tk. 11.15. Wx 



I. INTRODUCTION 



A variety of novel features has spurred interest in low-dimensional QED for many decades[Tj]. For example, high 
temperature QCD can be represented as the dimensionally reduced QCD3. If the number of quark flavours (Nf) is 
large, the nonabelian behaviour of the theory is suppressed and it may be approximated as quantum electrodynamics 
' in three dimensions (QED3)@. Massless QED3 in the large Nf limit generates dynamical fermion masses that are 
suppressed exponentially in the fermion number. Thus this theory illustrates how large mass hierarchies can be 
O : dynamically generated [3j, which is of interest to BSM physics. 

More recently, QED3 has been used as a model field theory for three dimensional condensed matter systems. 
Examples include applications to high T c superconductors, where the relevant dynamics is thought to be isolated to 
copper-oxygen planes in the cuprateQ. It is also considered as a gauge formulation of the 2+1 dimensional Heisenberg 
spin modelQ, a possible model for graphene[6], and quantum versions of spin-iceQ. 

When coupled to Nf massless four-component fermions, QED3 exhibits a U(2Nf) symmetry that can be broken 
to U(Nf) x U(Nf), which is the three dimensional analogue of chiral symmetry breaking. It is widely agreed that 
symmetry breaking occurs for low Nf. More interesting is the behaviour as the number of fermions becomes large 
| where it is possible that the interaction becomes sufficiently screened that vacuum condensation no longer occurs. 
Indeed, Appelquist et al. have used the large Nf limit with additional approximations to argue that there is a 
critical number of flavours, N+ rs 3.5 above which the theory remains in the symmetric phase. Furthermore, the 
lack of massless scalar bound states at the critical coupling has been used to argue that the chiral restoration phase 
transition is not second order, but is of a novel type[l]. 

The belief that a critical number of fermion flavours exists is not without controversy. For example, Pennington 
and Walsh assert that wavefunction renormalisation (which is null as Nf — > oo) is central to determining the details 
of symmetry breaking. Their solution of a truncation of the Schwinger-Dyson equations reveals that, in fact, chiral 
; I ■ symmetry is never restored[9(. In addition, Pisarski has used the renormalisation group to argue that fermion mass 
generation occurs for all values of [l(J. In spite of these claims, more recent Schwinger-Dyson computations again 
find a critical value of N f with iV+ rj 3.5 [H Ell- 
in principle lattice gauge theory can determine the nonperturbative properties of QED3. In practice, lattice 
computations have been hampered by the large ratio of the coupling and dynamical mass scales. Compounding the 
difficulty is the long range nature of the force, which leads to large finite volume effects. Nevertheless, computations 
with large lattices have been made and find N* ~ 1.5 [13, EH- It appears that the discrepancy is due to the necessity 
of extremely large lattices to achieve the continuum limit. This has been explicitly demonstrated by Schwinger-Dyson 
calculations with QED3 on a torus, which also show that quantities can appear to falsely converge at smaller lattice 
sizes [Hi. 

A novel feature of QED3 is that it is possible to introduce a topological Chern-Simons-like photon mass term to 
the theory 16]. This term breaks parity and time reversal symmetries. It is also possible to formulate the theory with 
two-component fermions. In this case a nonzero photon mass induces a finite fermion mass at one-loop (and vice 
versa) [H El- This raises the interesting possibility that parity symmetry can be spontaneously broken in the massless 
theory. This question was first examined by Appelquist et al. many years ago[18|. They concluded that finite fermion 
masses were dynamically generated, but that these masses appear in pairs of opposite sign, thereby maintaining 
the parity symmetry of the vacuum and a massless photon, in agreement with a general argument of Vafa and 
Witten[l9(. Their conclusions were based on an analytic examination of the Schwinger-Dyson equations in rainbow- 
ladder approximation. In addition the authors assumed the large Nf limit, that there is no fermion wavefunction 
renormalisation, that the fermion self energy is constant at low momentum, and that it is valid to truncate integrals 
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at the scale e 2 . Their general conclusions were confirmed by other groups with similar methods j20|. Nevertheless, a 
surprising counter claim exists, namely Hoshino and Matsuvma(2l| obtained a nontrivial parity violating numerical 
solution to the Schwinger-Dyson equations for Nf = 1. This analysis employed the rainbow-ladder approximation 
and a bare photon propagator. There appears to be, furthermore, a community that believes that spontaneous parity 
symmetry breaking is viable 22j. 

We revisit the massless two-component theory in an attempt to resolve the discrepancy mentioned above and to 
examine the robustness of the parity symmetry preserving scenario under a variety of truncation schemes. Since 
the two-component theory in a reflection invariant vacuum is equivalent to the four-component theory with Nf / 2 
fermions, we can also study the possible existence and stability of the chiral symmetry restoration transition. Of 
course one generally seeks stability with respect to truncation variation as a sign of robustness. Unfortunately, this 
can be difficult to assess in a rather small truncation model space. Gauge invariance is also a useful metric if suitable 
observables can be found. The situation is alleviated somewhat by a theorem due to Coleman and Hill, which gives 
the exact result for the parity- violating portion of the photon self energy j23|. We will use this to check the efficacy of 
various truncations in the following. 

Our investigation will utilise the Schwinger-Dyson equations truncated at the two-point function level with various 
model fermion-photon vertices. Our most sophisticated model will employ the Ball-Chiu longitudinal vertex with 
the Curtis-Pennington transverse vertex, putting the computation on a par with the most recent investigations of 
four-component QED3jl2|. We will confirm that nontrivial fermion masses are generated in opposite sign pairs, and 
therefore the vacuum remains reflection invariant. Furthermore, we confirm that chiral symmetry is spontaneously 
broken at low Nf and is restored at a critical number of flavours that depends weakly on the vertex Ansatz. Finally 
we establish that the Coleman-Hill theorem is satisfied reasonably well and that this agreement improves with the 
quality of the vertex, lending support to the conclusions presented here. 



We study massless abelian gauge theory in three dimensions with two-component fermions. The Lagrangian is 
taken to be 



The coupling e 2 has units of mass and the theory is superrenormalisable. As mentioned in the introduction, a fermion 
mass term generates a topological photon mass 



II. QED3 AND THE SCHWINGER-DYSON EQUATIONS 



C = -\f 2 + $(ip + em ~ ^(9 ■ A) 2 . 



(1) 



(2) 



via radiative corrections. Similarly, a nonzero value for fi will generate a fermion mass at one 
The full photon propagator is given by the expression 



loop mini. 




(3) 



where the projector is defined as 



(4) 



and the full photon self energy is parameterised as 



n M „0) = P 2 ^V n (P) + ie ^aP' 



n(p). 



(5) 



The precise statement of the Coleman-Hill theorem is 




(6) 
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to all orders in perturbation theory. This quantity is therefore interpreted as the topological photon mass. We have 
introduced a = e 2 /(Air). All subsequent numerical results are expressed in units of a. 

The proof of the theorem is relatively simple: if one considers the effective action obtained by integrating out 
fermions, gauge invariance implies that the leading behaviour of a general n-point function is 



r (n) (P!-fo) = o(pift). (7) 

Two external photon lines contribute 0(p 2 ) to any diagram beyond one-loop and so can not contribute to a function 
that is proportional to p. Nonzero contributions can only arise from graphs in which the two external photon lines 
end on a loop that has no other photon lines attached to it. Thus all corrections to the topological mass beyond one 
loop order vanish [30(] - It is also known that the two-loop correction to n(0) vanishes [25| . The two results together 
suggest that corrections to the topological mass vanish order-by-order in perturbation theory. We will shortly develop 
expressions for II and II in the Schwinger-Dyson formalism. Since an arbitrary truncation of the Schwinger-Dyson 
equations sums diagrams in a way that need not be consistent with perturbation theory or gauge invariance, one 
cannot expect that the Coleman-Hill constraint on 11(0) holds. Deviations from Eq. [5] can then serve as a useful 
diagnostic for the accuracy of a given truncation scheme. 




FIG. 1: Schwinger-Dyson Equations. Solid circles represent full propagators. The open circles represent a model vertex. 



The Schwinger-Dyson equations for the two-point functions are shown in Fig. [T] Solutions to these equations yield 
the full photon propagator, parameterised in Eq. [3] and the full fermion propagator, defined by 

The Schwinger-Dyson equations have been truncated by assuming a model form for the vertex (denoted with a small 
open circle). Typical model vertices include the rainbow ladder approximation: 

iT^ RL (k,p)=^, (9) 

the central Ball-Chiu vertex 

iT» CBC {k,p) = \{A{k) + A{p)) 1 », (10) 

or the Ball-Chiu vertex 



iY» BC {k, P ) = \ (aw + a( P )) r + l ^^h t+tw+pn ^L^ y+z). (ii) 

The Ball-Chiu vertex is the unique form of the longitudinal portion of the vertex that is consistent with the Ward- 
Takahashi identity and is free of kinematic singularities (2(|. Use of the Ball-Chiu vertex is a necessary, but not 
sufficient, condition for gauge invariance of the solution. 

The transverse portion of the full vertex remains unspecified. Curtis and Pennington have used multiplicative 



renormalisability to argue that the most important transverse term is[27| 



«t&p(M = \ m d{k ^ ] [r(k 2 - p 2 ) - (k + pT (je - i>)] 



(12) 



with 



d(k,p) 



2\2 



(fc 2 - p 2 ) 2 + (M(fc) 2 + M{p) 2 ) 
k 2 + p 2 



(13) 
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where M — B/A is the mass term of the full propagator. We shall refer to these truncations as RL (rainbow-ladder), 
CBC (central Ball-Chiu), BC (Ball-Chiu), or CP (following Ref. Q3, this is the Ball-Chiu plus Curtis-Pennington 
vertex in the fermion propagator and the Ball-Chiu vertex in the photon propagator). 
In rainbow-ladder approximation the Schwinger-Dyson equations for the ith fermion are 



' (2tt)3 { (A?q*-B?)D(K) + K 2 (A 2 q 2 — Bf)D(K) 2K 2 (Afq 2 - Bf) J ^ 



and 



2 .p-K(n-U), 1 



We have defined K = p — q and 



D(K) = K 2 (l - U(K)) 2 - (n - fl(K)) 2 (16) 

in these expressions. 

The rainbow-ladder expressions for the scalar photon functions are 



2 m 2, 9 , 2^ f d 3 q Mq 2 )MQ 2 ) ( o q-pQ-p \ f . 

P U(p ) = 2*e g J (AV _ B 2 )( M Q) 2 Q 2 - Bi(Q) 2 ) V 3 ^^J (17) 

and 

2tt ^_ o-2^ /• rf 3 g g-P^(g 2 )^(g 2 )- g -M,(g 2 )^(Q 2 ) 

pi^pj Ae ^y(27r)3 (AV-^X^fOT-BiW) 3 ) ' 1 j 

where Q = p + q. The expressions with the BC or CP vertices are more complicated and not very illuminating and 
hence are given in the Appendix. 

Note that the anomalous photon function n arises in the equations for At and Bi, thus hidden parity symmetry 
can affect chiral symmetry breaking and vice versa. Furthermore, (for fi = 0) the relation n[— B] — —H[B] implies 
that the Schwinger-Dyson equations are invariant under B — > —B. This remains true with the Ball-Chiu and Curtis- 
Pennington vertices. 

Finally, computing perturbatively with dimensional regularisation yields the following results for the one-loop scalar 
photon functions: 



n = -2a V / 



\Jrni — x(l — x)p z 



and 

iff 

(20) 



"f pi 
II = a > / dx 



\/m 2 — x(l — x)p 2 



Thus n(0) = —a/ (3m)Nf and n(0) is given by Eq. |6] We have confirmed that these results are recovered numerically 
when A = 1 and B — m. 

We shall consider two order parameters in the following, the dynamical mass 



M(p) = B(p)/A(p) 



(21) 



5 

and the fermion condensate (ijjip). Direct computation yields the Euclidean result 



Alternatively, examination of Eq. [14] reveals 

B(p 2 ) -> 2ttq^-±^ W)- (23) 

We stress that this result is true regardless of the vertex approximations made. Its general validity derives from the 
operator product expansion [28|. 

As mentioned in the introduction, Appelquist et al. have analysed the Schwinger-Dyson equations in the large Nf 
limit and obtained a critical number of fermion flavours, N+ rs = 64/ V 2 for two-component fermions, above which 
chiral symmetry is restored. Of course, one can criticise this computation for deri ving a relatively small value of A* 
in the large Nf limit. More recent versions of this calculation are made on Refs. and [l2|]. The result is of the 
form 



B(0) oc N f exp f 2n | . (24) 

As mentioned, this result is not entirely accepted; we shall analyse its accuracy numerically in the next section. 

Finally, Pennington and Walsh have used renormalisation group arguments to suggest that the wavefunction renor- 
malisation behaves asQ 



p 2 _ l _ lrn 2\*/**'*tf 



MP) ~ ) • (25) 

We consider the case m = so this equation might imply that A runs to zero with momentum. However the dynamical 
fermion mass is generically expected to mimic the bare fermion mass, and one can anticipate that A is finite at the 
origin but approaches zero as Nf is increased past A*. 

Finally, we choose to work in Landau gauge, which is widely regarded as providing the most reliable results when 
working at this level of truncation. For detailed discussions of gauge variance we refer the reader to Ref. [29| . 



III. NUMERICS AND RESULTS 



A. Numerical Methods 



The generic numerical problem we seek to solve is of the form X = F[X] where X represents the unknown functions 
A and B and F is a functional denoting the nonlinear integrals on the right hand side of Eqs. 1 141 151 Notice that the 
photon functions are relationships, not equations to be solved. We employed a number of strategies to solve these 
equations, including: 

(i) natural iteration. 



= F[JTW] (26) 

(ii) implicit iteration 



X {n) = F[X (n+1) ] (27) 

(iii) iteration with over-relaxation 



= (1 - + uF[xW] 



(28) 
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(iv) Newton-Raphson iteration 



X(«+!) = (l 



SF_ 

sx 



[xW]) _1 (F[xW]-jrW) 



(29) 



(v) minimisation 



min||X-F[X]||. 



(30) 



In addition, a variety of representations of the functions A and B were attempted including the point basis, modified 
Legendre polynomials, and a quadratic+inverse quadratic spline. We found that natural iteration worked quite well 
for all vertices except Ball-Chiu. Of course simply iterating need not yield a solution at all, rather the algorithm 
could approach a limit cycle, or could even go through a period doubling transition to a chaotic regime. Newton- 
Raphson iteration is both faster and more stable (but not always) than natural iteration. Over-relaxation works quite 
well if Gauss-Seidel updates are used. As with elliptic partial differential equations, one must tune the value of uj 
carefully. We did not find cases where implicit iteration was useful. Minimisation is always stable, but can be very 
slow, especially when a large number of degrees of freedom are being varied. Unfortunately, the method can easily 
converge to one of a great many local minima. Worse, these minima need not be close to an actual solution to the 
Schwinger-Dyson equations and care is required. 



We seek to confirm the reflection symmetry breaking pattern suggested by Appelquist et ai. Thus if parity symmetry 
is maintained one expects either no chiral symmetry breaking or chiral symmetry breaking with Nf/2 fermions of 
mass M and Nf/2 fermions of mass —M. If reflection symmetry is broken then one expects Nf fermions with mass 
M (or — M). (Other patterns are possible; we restrict attention to these three). 

We have performed computations assuming a reflection symmetry breaking pattern of (M ,M) with all vertex 
models. The natural iteration algorithm approaches a limit cycle of the type (A,B) —> (A,—B). While this could 
prove to be a nontrivial solution (due to the B reflection symmetry discussed above), the error \\X — F[X]\ \ does not 
decrease with iteration. This suggests that the limit cycle does not represent an actual solution to the Schwinger- 
Dyson equations. Of course, the lack of convergence of a natural iteration algorithm does not prove that a solution 
does not exist. We have therefore repeated the computation with the minimisation algorithm. In this case A evolves 
to a smooth function while B approaches zero (with substantial noise). Similar behaviour was seen in all vertex 
models. Thus it appears that a nontrivial reflection-breaking solution to the Schwinger-Dyson equations does not 
exist, as suggested in Refs.[l8j], [20T |. This computation is the first extension of those conclusions to the full two-point 
function Schwinger-Dyson equations with vertex models of varying sophistication. 



Recall that Ref. [2l| found a parity-violating solution. Their computation was made in the quenched rainbow- 
ladder approximation to the Schwinger-Dyson equations. We, in fact, confirm their result; however, it is clear that 
the quenched approximation is unsuitable for investigating properties of QED3. The photon propagator is changed 
drastically in the infrared, even in perturbation theory, and this must be accounted for in the formalism. 



We have argued that the value of 11(0) can be used to test the accuracy of truncation schemes in the Schwinger- 
Dyson equations because its value is fixed by Eq. [5] to all orders. We display 11(0) as a function of the number of 
fermions for the RL, CBC, BC, and CP Ansatze in Fig. [2] As hoped, more sophisticated vertex Ansatze maintain 
the Coleman-Hill result more accurately and over a larger range of Nf. The CP vertex remains within 20% of unity 
over the range of Nf considered. This is a highly nontrivial phenomenon: as Nf increases B becomes exponentially 
small, driving II to also be exponentially small (see Eq. I18p . This is countered by an increasingly exposed infrared 
divergence in the integral. It is remarkable that the two effects cancel so precisely when the CBC, BC, or CP Ansatze 
are employed. Because of this, we estimate that the truncation error in our CP results is less than 20% over the range 
of Nf considered here. 



B. Results 



1. Parity Symmetry 




2. The Coleman-Htll Constraint 
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1 2 3 4 5 6 



Nf 

FIG. 2: fl(0) vs. Nf for various models. From top to bottom at Nf = 5 these are: RL, CBC, CP, BC. 



3. Chiral Symmetry 

As expected, we find chiral symmetry breaking solutions for small Nf for all vertex models. The A/ -behaviour of 
the solutions is explored in Figs. [3]and|4j One sees that B(p) drops precipitously with increasing number of fermion 
flavours. Whether B vanishes at a finite value of Nf will be discussed below. First we note that the value of A drops 
near the origin and may be approaching a limit function with A(0) — 0. Does this happen at a finite Nf? A similar 
behaviour is seen for II, again, with the possibility that an infrared singular limit function is approached at a finite 
value of Nf. These figures have been generated with the CBC vertex, qualitatively, and in most parts quantitatively, 
similar results arc obtained with the other truncations. 
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FIG. 3: A vs. p for various Nf (L). B vs. p for various Nf (R). These curves are obtained in the CBC truncation and correspond 
to Nf — 1, 2, 3, 4, 5, and 5.75 from top to bottom. 



We examine the putative phase transition more closely by plotting the fermion condensate versus Nf (Fig. [5J. 
We have confirmed that the directly measured condensate agrees with that obtained from Eq. [521 to high accuracy. 
Results for the RL, CBC, BC, and CP models are shown in Fig. [5l Also shown are fits with the Ansatz 



{ ^ m) = aN ^{ ^N~/N f -l ) 



(31) 
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FIG. 4: Vacuum Polarisation vs. momentum for various Nf. These curves are obtained in the CBC truncation and correspond 
to Nf = 1, 2, 3, 4, 5, and 5.75 from bottom to top. 
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FIG. 5: The condensate for (top to bottom) BC, CP, RL, and CBC models, with fit functions (described in the text). 

It is quite difficult to obtain these results for large values of Nf. extrapolation in the number of integration points 
is necessary for Nf > 4. Our most accurate results are for the CBC vertex and are in excellent agreement with the 
Ansatz of Eq. 13 X I implying strongly that a chiral restoration phase transition does indeed occur. An Ansatz like 
00 ex P( — Nf) also does a reasonable job fitting the CBC points, and of course does not exhibit a phase transition. 
Nevertheless the quality of this fit is significantly lower than that of Eq. |3"T1 

Remarkably, a simple adjustment of A* in Eq. [31] provides excellent fits to the RL, BC, and CP data as well. We 
obtain a = 1.8 and 



A* (CBC) = 1.00 ■ Nf 



(32) 



A*(RL) = 1.10 • Nf 



(33) 



A*(CP) = A*(BC) = 1.21 ■ Nf. 



(34) 



Although we can be less certain that the BC and CP models exhibit symmetry restoration, the good fits, the similarity 
of all model results, and the theoretical motivation all make this very likely in our opinion. Our numerical results 
appear to be in agreement with earlier results [ill . [l2| . 
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IV. CONCLUSIONS 



Previous arguments for the lack of reflection symmetry breaking in two-component QED3 relied on the large Nf 
expansion and several simplifying assumptions. We have extended the analysis to finite Nf with several model vertices, 
including the Ball-Chiu-Curtis-Pennington model, which is the most sophisticated model of QED3 to date. We cannot 
find nontrivial solutions for the fermion propagator with several numerical techniques, making it very likely that parity 
symmetry is maintained in two-component QED3. 

A novel feature of this analysis has been the use of the Coleman- Hill theorem as a diagnostic for truncation accuracy. 
General arguments imply that truncations that obey the Ward-Takahashi identity should perform better than those 
that do not. We have confirmed this. We observe that deviations from the exact value of 11(0) indicate truncation 
errors of approximately 20% over a large range of Nf, so there is some hope that the conclusions obtained here are of 
reasonable fidelity. 

Parity-symmetric two-component QED3 is equivalent to four-component QED3, permitting us to examine the is- 
sue of chiral symmetry breaking and restoration. The difference between symmetry restoration and exponentially 
suppressed symmetry breaking is necessarily small, and numerically distinguishing these scenarios can be very chal- 
lenging. While not definitive, our results agree with earlier computations in the Schwinger-Dyson formalism and 
hence support the original simple arguments of Appelquist et ai. Furthermore, it is likely that the transition occurs 
in all the model truncations considered here. The value of the critical number of fermions varies by model, but all 
are surprisingly close to the original estimate. 



Acknowledgments 



This research was supported by the U.S. Department of Energy under contract DE-FG02-00ER41 f 35 and an Andrew 
W. Mellon Predoctoral Fellowship. 



V. APPENDIX 

The truncated Schwinger-Dyson equations for the fermion propagator with the BC+CP vertex arc as follows. 



o / d 3 q 

B(p)=m- le I 



A + d2 AA AB . d2 



04 



(35) 



and 



At ^ 1 f d3 1 

A(p) = l + i- / — ^ 

P J (2tt) 



A + d2 AA AB 



d2 



-p-a 1 + —a 2 + —a 3l 



0:4 



(36) 



Subexpressions are 



A=-(A(p)+A(q)) 



AA 



1 A(q) - A(p) 

2 q 2 — p 2 



(37) 
(38) 



AB 



B(q) ~ B(p) 
q 2 — p 2 



D = A( g y q 2 - B(q) 



d2(p, q) 



l (A(q)-A(p))(q 2 -p 2 ) 
2 d{q,p) 



(39) 
(40) 
(41) 



See Eq. EH for d(q,p). 
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The coefficients are 



Pi = 2aB - 2Abq ■ K + £B/K 2 



p 2 = aB(Q 2 -Q-KQ-K) + £B(Q ■ K) 2 / K 4 - Ab{p ■ Qq ■ K - p ■ Kq ■ Q) 
p 3 = Aa(Q ■ q - Q ■ Kq ■ K) + £AQ ■ Kq ■ K/K 4 
p 4 = -£B(p 2 - q 2 )/K 2 - 2Ab(p 2 q 2 - (p ■ qf ) 
on = -2Aap ■ Kq ■ K + £,A(2p ■ Kq ■ K - p ■ qK 2 ) / K 4 + 2Bbp ■ K 

a 2 = ^Aa(q 2 + p 2 )(Q 2 - (Q ■ K) 2 ) + £A(Q ■ K) 2 p ■ q/K 4 - Bb{p ■ Qq ■ K - p ■ Kq ■ Q) 

a 3 = -Ab(p -Qq- K -p- Kq-Q) + aB(p Q-p- KQ ■ K) + £Bp ■ KQ ■ K/K 4 



a 4 = -2Bb(p 2 q 2 -{p-q) 2 ) + A 

+ q-K(p-Q-p-K(p 2 -q 2 )/K 2 )}- 

i(2p ■ Kq ■ K(p 2 - q 2 ) - K 2 p ■ q(p 2 - q 2 )) / K 4 

1 - U(K) 



a[p-K{q-Q-q- K(p 2 - q 2 )/K 2 ) + 



K 2 (1-U{K)) 2 -T1(K) 2 



fl(K) 



K 2 (1-U(K)) 



K»=p^- g" Q» = p» + <f 
Ball-Chiu expressions for the photon scalar functions are 



and 



p 2 Tl(p) = 2ie< 



p 2 fl(p) = -2ie 2 



d?q 
(2irf 



d?q 
(2tt) 3 



A AA AB 
D* 1 + ~D n2 + -W* 3 



A AA AB 



In this case 

D = [A(qfq 2 - B(q) 2 ] [A(Q) 2 Q 2 - B(Q) 2 } 
and the subexpressions are given by (note that the arguments of A, AA and AB are now Q and q) 

Tri=A(q)A(Q)(q-Q-3q-pQ-p) 

2tt 2 = A(q)A(Q) [2q ■ (q + Q)Q ■ (Q + q) - (Q + qfq ■ Q - 3q ■ p(Q + q) ■ pQ ■ (Q + q) + 
+ 3(Q + q) ■ p(Q + q) - pq-Q - 3Q ■ p(Q + q) ■ pq ■ (Q + q)] + 
+ B(q)B(Q) [(Q + q) 2 - 3{Q + q) ■ p(Q + q) ■ p] 



(42) 
(43) 
(44) 
(45) 
(46) 

(47) 

(48) 



(49) 
(50) 

(51) 
(52) 

(53) 

(54) 

(55) 
(56) 

(57) 
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2n 3 = A(q)B(Q)(q .( q + Q)-3q- p(Q + q)-p) 



+ 



B(q)A(Q)(Q ■(Q + q)-3Q-p(Q + q)-p 



(58) 



fri = B{q)A(Q)p ■ Q 



A(q)B{Q)p ■ q 



(59) 



n 2 = (A(q)B(Q) + A(Q)B(q)) ■ {q 2 p 2 - (q ■ pf) 



(60) 



2n 3 = A(q)A(Q)(q ■ (Q + q)p 2 - p ■ (Q + q)p ■ q) 



(61) 
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